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induced by heavy scalar field in adjoint representation at N = oo for arbitrary dimen- 
sion D. The nonlinear integral equation for the gauge invariant density of eigenvalues 
of the vacuum average of the scalar field is derived. In the continuum limit, the density 



We find the exact solution of a recently proposed model of the lattice gauge theory 
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1 Introduction 

Recently, the new lattice gauge model was suggested by Kazakov and myself Q. We observed, 
that the Yang-Mills term in the Action could be induced by massive scalar field in adjoint 
representation, in the same spirit as actual QCD is induced by some yet unknown quark 
and gluon constituents. The origin and properties of these constituents is irrelevant for the 
purposes of QCD, in virtue of renormalizability. 

The basic property of induced gauge models is the scaling law for the physical mass scale 
M as a function of the constituent mass mo 

M oc (m — m c ) 7 (1) 

where m c is some critical value, depending upon the regularization scheme, and the critical 
index 7 depends upon the universal characteristics of the constituent field, such as spin and 
number of components. 

This scaling law comes about when one follows the running gauge coupling constant 
P( r ) = ~4ns from the lattice spacing scales r = r where we assume (3{tq) = all the way 
to the hadronic distances r 2 ~ jj , where again /3(r 2 ) ~ due to confinement. Along the 
way, /3(r) first goes up as c\ In — due to screening by the constituent field, until we reach the 
scales r\ such that ^(wq - m^) ~ 1 where the screening turns off, leaving us with some large 
constant @{r\) = Ciln^ 1 . This constant serves as the bare coupling for the gauge theory, 
with effective ultraviolet cutoff at r\. The rest of the story is well known: the /3(r) goes down 
as 0(ri) — c 2 In — due to famous antiscreening from induced selfinteraction of the gauge field, 
all the way to the physical scale r 2 . Comparing two parts of this A shaped curve in the lnr, /3 
plane we find c\ In — = c 2 In — which yields the scaling law ([!]). 

The above arguments provide only a rough estimate of induced coupling. There are 
feedback effects from the hard gluons, as well, as some induced selfinteraction of the inducing 
field. If we go back to the cutoff scales from the scale r\ with some unknown constant f3\ we 
would have to compute the renormalization of the effective </> 4 theory coupled to gauge field. 
The peculiarity of this theory is that the scalar particle is very heavy, not like the ordinary 
Coleman- Weinberg model. Clearly, the effective 4 coupling cannot be small, otherwise 
the screening of charge from scalar particle would not overcome the antiscreening from the 
gluons. Indeed, as we shall see in that last Section, it must be adjusted to some critical 
value, in order to get the scaling law. 

We already know the example of the heavy particle theory which in a large iV limit 
becomes the string theory rather than particle theory. This is the same scalar field in 
adjoint representation of SU(N) but without the gauge field. If we take the large N limit, 
then we are left with planar graphs, which have singularity at some value of the coupling 
constant. So, if we adjust the scalar coupling to approach this large N singularity, instead 
of adjusting the bare mass to get the small physical mass, then the planar graphs are known 
to condense to a fishnet with lattice scale cells, leading us to the string theory. 

The problem with the latter model is well known, it does not seem to have a continuum 
limit at D > 1. The model we are considering now, is slightly different. If you look at the 
strong coupling expansion of the induced gauge theory, by integrating the scalars first, and 



representing the effective gauge action as sum over the scalar loops |1|], you find the sum over 
closed surfaces made of these loops glued together by sides. The typical size of the loop is 
ri, in above notations. In the large N limit we expect only the planar surfaces to survive, 
and produce some kind of string theory, with possible extra internal degrees of freedom at 
the world sheet. 

The resulting theory must be universal, as we expect it to fall into one of the fixed 
points. In four dimensions there is only one nontrivial fixed point known, the QCD. So, the 
crucial question is whether the solution of particular induced gauge model is nontrivial. This 
question is a dynamical one, and could be answered only by exact solution. 

Our best bet here is the scalar model in adjoint representation. The unique property of 
this model is the possibility to diagonalize the scalar field by the gauge transformation 

$(x) = n(x)<p(x)n\x); n\x)n(x) = 1 (2) 

so that the number of scalar degrees of freedom drops down from N 2 — 1 components of 
the traceless Hermitean matrix $(x) to its iV — 1 eigenvalues <f)(x) = diag(<j) a (x)),a = 
1, . . . N, 2 <Pa — 0. These gauge invariant eigenvalues serve as a master field: they freeze at 
a certain vacuum average, with fluctuations Sep ~ -k- 

The active degrees of freedom in the large N limit are given by the gauge fields U^x), 
associated with the links of the lattice. The analysis of |l[] shows that the usual scalar particle 
loops in external gauge field within the strong coupling ( random walk) expansion can be 
reproduced by the U(N) group integrations over the gauge field for the frozen master field. 
The summations over eigenvalues, distributed in the strong coupling expansion according to 
the Wigner semi circle law 

p{4>) oc y/a 2 - 4> 2 (3) 

yields in the large N limit the same results as Gaussian integrations over matrix field $(x). 
If we accept this scenario, then the model can be dramatically simplified by integration 
over gauge fields independently at each link. The catch is that this integral || 

/(*,*) = (DUaviNtcSmtf) oc ^i^Ml (4 ) 

where 

a(0) = no* - <f>j) ( 5 ) 

is the Vandermonde determinant, becomes at iV — > oo nontrivial statistical system in external 
field, given by eigenvalue distribution 

p*(v) = ^E^-^i p*(y) = j^Y. 5 ^-^) ( 6 ) 

■'' a ■'"a 

This external field problem should be solved together with the saddle point equation 

"ain/($,^) 



(p a — jJL I 7V->oo j\ 



(7) 



where D = 4 is dimension of space, which we shall keep arbitrary in the rest of this paper. 
We introduced the general scalar particle potential 

U{4>) = \ml<f + iA0 4 + . . . (8) 

rather than the pure mass term considered in [p]]. The higher terms would be needed to get 
to the correct fixed point of our model. Note, that we defined the mass term with traditional 
| factor, which was skipped in Jl|. 

In the previous paper |l.| existence of two phases of this system was conjectured. The 
so called weak coupling phase, with a gap in the eigenvalue distribution reduces to the two 
matrix model 

/($,*) oc Jd N xd N yA(x)A(y)exp U^^Vi ~ ^i(^) ~ ^(w)J (9) 

with potentials 

V$(x) = / d/ip$( / u)ln(/i — x); Vq(y) = / dup^(u) \n(u — y) (10) 



which model was solved exactly ai N = oo. 

In the strong coupling phase, where there is no gap, the same representation @ is 
formally valid, but potentials are singular, which leads to drastic changes in the large N 
limit. The index sums cannot be replaced by the integrals in equations of the orthogonal 
polynomial method, but the individual terms become singular and produce the necessary 
factor of N. We managed to generate several terms of strong coupling expansion from the 
two matrix model, but could not find the exact solution. 

This is a pity, since the induced QCD scenario corresponds to the strong coupling phase, 
where the distribution of eigenvalues is qualitatively the same as the semi-circle law. There 
are also some internal inconsistencies in the weak coupling solution]!]], in particular, there 
are no physical excitations in the vacuum, as the hopping term in the effective Lagrangean 
for the vacuum fluctuations vanishes in the large N limit. 

In the present paper we find the relevant strong coupling solution by completely a different 
method. The solution of the lattice model is described in Section 2, where we derive the 
lattice version of the equation for the density of eigenvalues of the scalar field. The critical 
phenomena are investigated in Section 3, where we find a remarkably simple local version of 
this equation and solve it at the critical point. 

2 Exact Solution in the Large N Limit 

We are going to apply to the one link integral the good old Schwinger- Dyson equations. The 
traces of properly normalized derivatives with respect to $ matrix 

< v *'« s juk, (11 > 
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act as traces of \l/ matrix on our integral 

tr (Vg-^ n )/($,^) = (12) 

as it follows from the unitarity condition. Let us introduce the matrix function 

F($) = V* In /($,#) (13) 

and rewrite above equation as follows 

tr ((F($) + V$) n * 1 - ^ n ) = (14) 

Now, summing the geometric series, we define the matrix 

G A ($) = (A - F($) - V^r 1 * 1 (15) 

which satisfies the linear differential equation 

V„G A ($) = -1 + (A - F($))G A ($) (16) 

and extra condition 

tr (<*<*) + wh) = ° (17) 

The differential equation allows us to express G in terms of F at fixed spectral parameter 
A, after which the last condition would provide equation for F. One may readily check 
that these equations work in the strong coupling phase, where the function F(<f>) expands in 
power series of its matrix argument. This expansion goes along with the expansion of G\(<&) 
in inverse powers of A, generating the moments relations (13). The expansion coefficients 
of F(<&) are obtained one after another from the moments relations. We checked that the 
same coefficients arise in the orthogonal polynomial treatment of the two matrix model in 
the strong coupling phase. The first several terms of this expansion were presented in refJT]. 
Now, let us take the large N limit. The key observation is that at large N one does not 
have to differentiate the traces, as that would not produce the necessary factor of N. Such 
factors are only produced by derivatives of the powers of the $ matrix 

n i 

fc=l N 

which terms remain finite in the large N limit, when this formula can be written in an 
integral form 

V*$" = J d/i/nM^- (19) 



This representation allows us to write the linear integral equation for function G\{y) 

I d W M GM ~ Gx{u) = -1 + (A - F{u))G x {v) (20) 

J a — v 



As for the extra condition (1171), it can be rewritten as follows 



JdfipMGM = Jdv^- = V>(\) (21) 

At the moment we are interested in the spatially homogeneous master field, so that \I/ = $, 
in which case we can omit the subscripts of the densities and potentials. 

We define these spectral integrals in the principal value sense, so that both F(/j,) and 
G\(fi) are real functions, when A , \x vary at the real axis. At A — > oo by construction 
G\(n) — > t- We shall assume that the eigenvalues are distributed along the whole real axis, 
this assumption will be justified later.[] 

The saddle point equation allows us to eliminate one of the unknown functions, by ex- 
pressing F{v) in terms of the potential 

F{u) = 2 l ^ ^ (22) 

When this expression is substituted back into fl2~0|), the terms combine nicely and we find 
the following singular integral equation 



1 + pf dfx P ^ GM = (A - R{is))G x {v); R{v) = F[y) + W{y) (23) 

J ii — v 

d as the boundary value problem for analytic function 

T x (z) = l+[dvM&M (24) 

J u — z 



At^pM 

This equation could be regarded as the boundary value problem for analytic function 

At)pM 

in the upper half of complex z plane, with usual analytic continuation to the whole plane by 

symmetry 

T x (z)=T x (z) (25) 

In virtue of symmetry of the eigenvalue distribution we could have mapped the upper half 
plane to the cut plane by conformal transformation £ = — z 2 but it appears to be simpler to 
work in the upper half plane, keeping in mind the symmetry of boundary values at the real 
axis. 

We shall assume the convergent integrals, in which case 

Tx(oo) = 1 (26) 

The next term of asymptotic expansion at z — > oo is also known in virtue of an extra 
condition (|2"1~D, namely 

TiW - 1 - ^ (27) 



lr This assumption, if fact will be valid only in the local limit, when the small vicinity of the origin in the 
eigenvalue distribution will be rescaled to the whole complex plane. In any particular lattice model, there 
might be finite support, but these effects are lattice artifacts, much in the same way, as the lattice momenta 
are distributed within the cell of inverse lattice. The eigenvalues have dimension of m^ D ~ l , which means 
that for D > 2 the cutoff for the eigenvalues grows as the power of the momentum cutoff. In the local limit 
both cutoffs disappear. 



The boundary problems of this kind can be solved in general. We express G in terms of 
discontinuity 

3T A (z/+) = np(ji)G x (n) (28) 

and find the following ratio of the complex conjugate boundary values of T at v± 

T\M A - R(u) + nxp(u) 



T\{yJ) A — R{y) — vnp(v 
The boundary problem has a unique solution 



_ , , f°° dv 

i\\z) = exp / — arctan 

\J-oo n(v — z) 



^ ^ (30) 



A - R(u) 



The ambiguity of the choice of the branch of arctangent is removed by an extra requirement 

that at A — > oo 

r lW - 1 - ^ (3D 



Comparing the first term at 2 — > oo with the asymptotics (27) we find the following 
integral equation for the potential 



-oo 



%y'(\)= I — arctan - ^7 x (32) 

j-oo 7r A — R[y) 

%V'(v + ) = -<k P {v) 

R{u) = ^U\u) + ^-WSV'iy) 

We shall refer to this equation as the Master Field Equation, or MFE. In Appendix we 
study general properties of this equation, assuming finite support. 

3 Critical Point 

Let us study the small vicinity of the origin in the MFE. The idea is, that the eigenvalues 
scale as positive power of mass, therefore, in the lattice units we are using, the continuum 
theory corresponds to infinitesimal vicinity of the origin. The normalization of the density is 
also the lattice artifact, as it follows from dimensional countuing. The internal nonlinearities 
in the scaling region, rather then the global normalization, would fix the scale of the density 
at the origin. 

With this picture in mind, let us restore the imaginary part of MFE from analyticity in 
the upper half plane 

r (:)= r±J 1 + ^M-) ,33) 

J-oo -in \ z — R[y) J 

In the lower half plane we would have — i instead of i in the same integral.^ 

2 Let us stress once again, that this formula is not valid for finite support of the eigenvalues. The correct 
formulas in this case are studied in Appendix, and this one can be obtained from the general formulas in the 
limit a — ► oo. 



To make it less surprizing, let us consider fi = — — as analytic function of £ = — ,2 2 , 
which it is, since the distribution of eigenvalues is symmetric. Then the upper half plane 
would map onto the first sheet of the Riemann surface, the real z axis being mapped onto the 
cut from — oo to in the ( plane, with negative (positive) z corresponding to upper(lower) 
side. 

Imaginary part of V'(z + ) is negative both for positive and negative z, and the real part 
changes sign with z, therefore the values of fi(C) a t the opposite sides of the cut are complex 
conjugate, as usual. The second sheet function, with i — > — i has the opposite sign of the 
imaginary part, so it matches the corresponding boundary values at the cut, as it should. 

We went into these boring details, because it took us some time to turn away from the 
usual picture with finite cut in the z plane and no singularity at infinity. We are going to 
use the old z variable in the rest of the paper, but restrict it to the upper half plane. 

Let us study the equation for the imaginary part — 7rp(A) of V'(A+). The imaginary part 

of MFE reads 

dv ( n 2 p 2 (v) 



ma)=/ Y^v + TBW=wy < 34 > 



oo 



Let us differentiate this equation, we find 



oo 



2 / \ 

p\\)=-pj dv— — P V — (35) 

(A - R[y)) ((A - R{v)) 2 + Ti l p\v)) 



oo 



Now, investigation of this equation reveals that critical phenomena take place when the 
renormalized scalar potential 

2«(0) = U{<f>) - Dtp 2 = l -m 2 <p 2 + i# 4 (36) 

goes to zero in the lattice units we are using, so that R{y) — > v.f\ 

The relevant eigenvalues would scale as some power of scalar mass v ~ m 7 , and density 
would scale as p ~ m 2+7 , the same as u'(<p) where the critical index 7 will be computed 
below. 

The scaling solution comes about in a tricky way. There are two relevant regions of 
integration in (^). In the region of \v — A| ~ m 7 , we could neglect p in the denominator, 
which yields the standard principle value integral. Note, that the size of the gap Si/ which is 
needed to define the principle value integral, could be any number in the interval m 2+7 ^C 
5v <^ m 7 , the result would not depend upon this gap then. 

There is another, much smaller region, that of \v — A| ~ m 7+2 ^C bv where we cannot 
neglect p but we can expand u'(u) and p[y) in Taylor series near v = A after which integration 
reduces to the residue in complex pole. The integral converges in this region, so that the 
boundary bv between regions drops from the answer. 



3 The remarkable fact, that in spite of all the strong fluctuations of the link variables, which were very 
far from continuum, we still get the same critical value rrv^, = 2D of the bare mass, as in asymptotically free 
gauge theory, would be discussed later. 



Collecting the terms, we find, that the term p'(X) cancels, so that we are left with the 
following integrodifferential equation 

/oo Cp'iu] 
dv — rr (37) 

-oo (A — V) 



r (A) = u'(X) + ^-^p r d\x ■ M 



D J-oo A — fj, 



/oo 2 (u) 
dfi (39) 

-oo Z — II 



The regular terms represent the contributions from the cutoff scales. As usual in the field 
theory, such contributions provide the contact terms, without singularities. Note, that this 
equation could be integrated twice in v at the expence of redefinition of the regular terms. 
On the left side we would obtain p(A)r(A), and on the right the denominator (A — z/) 3 would 
be replaced by 2 (A — v). 

This convenient method of elimination of regular terms is to introduce two analytic 
functions 

V(z) = y^u'(z) - V\z) (38) 

z — fl 
and note that at Qz — > +0 by construction 

%Q = (SV) 2 . (40) 

On the other hand, in virtue of the above equation for density 

^Q= l ~^^ (^ 2 ) (41) 

These two equations should be supplied with the symmetry conditions, namely, imaginary 
parts of these two functions at Qz — > +0 must be even functions of z, whereas real parts 
must be odd. The derivative of the real part of V(z) at the origin plays the role of effective 
scalar mass m 2 . 

This nonlinear boundary problem, which we call Local Master Field Equation, (LMFE) 
represents the main result of this work. At first glance it looks too simple to have nontrivial 
solution, but this is not true, as we found by analytical and numerical studies. 

We doubt that there could exist an exact solution in the general case, but precisely at 
the critical point, corresponding to the vanishing mass, such solution can be found in virtue 
of scale invariance. Let us take the powerlike Anzatz 

V(z) = V z°; Q(z) = q z 2a ; a = 1 + - (42) 

7 

The symmetry comditions would be satisfied provided 

p = ipo exp ( — lira J ; q = iq exp (— ina) (43) 



where po and q$ are real numbers. Now, the powers of z would cancel in LMFE, so that we 
may check it at z = 1 where 

9 V = po cos -vra; 3? 2 = pi sin Tea; SQ = g cos Tea; 3? Q = q sin ira (44) 

Comparing this with LMFE, we find 

D 2ix D 1 - D 2 

COS ™ = "3D = 2 ; C ° S 7 = 3D^2 ; * = ~^~ Po (45) 

The index 7(1?) takes special values 7(0) = 4,7(1) = 00,7(2) = 6 at low dimension, after 
which it smoothly goes to its asymptotic value 7(00) = 5 taking only transcendental values 
at integer dimensions. In particular, 7(4) = 5.41991, a(4) = 1.36901. Note, that unlike the 
singularity at D = 1 which was manifest in initial model, triviality of the theory at D = 2 
comes about in a very nontrivial way. In our opinion, this is an independent evidence for 
equivalence of this model to the gauge theory. 
The general scaling solution 

(46) 

starts from some constant at <f> = 0, then grows, approaching the above powerlike asymp- 
totics, p oc (p a at (j) ^> m 7 , so that the mass cancels. 

This is the magnified vicinity of the origin in complete spectrum of eigenvalues; ap- 
parently, the density reaches the maximum, and then starts decreasing outside the scaling 
region. The density at the origin goes to zero, so that at negative m? the gap would appear 
in the distribution (clearly, one cannot continue our irrational powers of m to negative m 2 
). In the forthcoming publication^] we solve this equation numerically at various D. 

As discussed in the previous paper M, the scaling laws in this model translate into loga- 
rithmic laws of original theory: for the effective induced gauge coupling of the form 




1 A 2 

-=■ ~ bN In ^5— 
gl ml - 2D 



— ~bN\n— — (47) 



the physical mass scale would be 




48irf_b 
11 



48tt 2 \ „ . (ml-2D\ 



n~u 




Should we identify this scale with the scale of our field {ml — 2D) 2 " 1 we would conclude 

that 

11 / 4tt \ 

"=4^ tr 1 ) (49 > 

but this seems to be too long a shot: one should investigate the physical mass spectrum 
before drawing any conclusions. 

9 
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A Investigation of the Master Field Equation 

Let us study the Master Filed Equation (|3~3"D assuming the finite support with the following 
Anzatz for the normalized spectral density 



OO 



1^1 — cos(2A;a) ,—, 

v = acos(a), < a < n, p(u) = - > £ k —— — , > 4 k = 1 (50) 

ttj^J osm(a) t^i 

This is just the Fourier expansion, with correct symmetry properties. The corresponding 
Fourier expansion for the real part of potential can be found by elementary integration 

sim/v raw f a a P( v ) ^ t sin W) / K1 N 

$lV{acos{(3)) = p / du — = }^Ck r— =r- 51 

J -a acos(p) — v ^i asm(p) 
and the potential as analytic function 



2k\ 



Substituting in (P3"|) the previous formula, and Fourier transforming, we find 



^/>W>°(^S) (^^T-^ (53) 

where we skipped for brevity the arguments of R, p. The contour encircles anticlockwise the 
interval —a, a. The last factor takes complex conjugate values exp(±2ik/3) at A = acos/3, 
which yields correct integrand after taking discontinuity. The term with p = is substructed, 
since this is how the branch of arctangent was defined. 
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The A integrand has only two singularities outside the contour, and decreases as integer 
negative power of A, so that there is no residue at infinity. The discontinuity along the cut 
from R — %p to R + ip is equal to n, and this integral is elementary. As a result we find 

(/ , \2fc-l / . \2fc+l' 
(c - vc 2 - i) (c - Vcp^T 
2,-T 1 A 2k +1 ' |-{g-»C) < 54 > 

where 

C = —^-L Ll—L = -JL cos a H — — S^V^' (a cos a) SsV' (a cos a + zO) (55) 

a 2D aD a 

We restricted ourselves here to pure mass term, without interaction terms in the scalar 
potential. 

Summing up the equation for £& from k — 1 to oo , using normalization condition ([50]) 
on the left, and summing up the series on the right, we find the following equation for a 

4= [* da sma% (c - VC 2 ^! + /(KC) 2 -l| (56) 

But the same normalization condition yields the simpler equation 

TV f n 

— - = / da sin a S C (57) 

a 2 Jo 

which means that contribution from the square roots should cancel. This relation can be 
written as follows 

du 1 — tJ" 



The explicit expression for the analytic potential ([52]) follows from above representation 



2 „2 



V\z) = \U- v^ 2 "^ 2 ) - aU ( Z — U " ) (59) 

Ll \ (Jj J 

where 

VS; io 2ttz V V ^ (C + VC 2 - 1-0/ 

Now, using (]5^) we can write down the integral equation for C(a). 

_, , a 2 m 2 -4 a 2 mg + 4D-2 2D-1,./1\ 1,,, x 

which should be solved together with (|58|) . 
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The semi circle law at ttiq — > oo corresponds to the following 

a 2 - 4- C(a) -> -^ (62) 

Wq a; 

since W(C) -> as C~ 3 at large C. The function C{ui) becomes analytic at infinity, which 
allows us to expand the contour and in (]58f) to infinity, where expansion goes in even inverse 
powers of u), so that there is no residue. As for the integral with the real part, here we cannot 
expand the contour, but the integral is small by itself, as the imaginary part of this square 
root does not vanish only in the small region of angles near |-7r. 

We solved these equations at D = 4 numerically by iterations, taking as unknowns the 
values of C at M = 4096 of points in the first quadrant 

C m = C [a = 2M + 1 ) ; < m < M; C 2M+ i- m = -C m (63) 

The integral was approximated by sum over these equidistant points, which corresponds to 
the Gauss-Tchebyshev method. 

As expected, at large ttt-q we reproduced the semi circle, moreover, it was difficult to get 
anything else. The nontrivial term U in our equation was practically negligible, unless we 
adjusted the bare mass so, that the derivative $IV"(0) vanished. Then we started seeing 
interesting phenomena, but in order to compare them with our local equations, much better 
resolution would be needed. According to our analysis in the Section 3, there are two different 
small scales involved, and there are strong cancellations. In such situation, the approximation 
of the integral by sum over equidistant points is no longer valid. The corresponding numerical 
solution is studied in the forthcoming workpfl 
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